High fidelity quantum gates via dynamical decoupling 
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Realizing the theoretical promise of quantum computers will require overcoming decoherence. Here we 
demonstrate numerically that high fidelity quantum gates are possible within a framework of quantum dynami- 
cal decoupling. Orders of magnitude improvement in the fidelities of a universal set of quantum gates, relative 
to unprotected evolution, is achieved over a broad range of system-environment coupling strengths, using recur- 
sively constructed (concatenated) dynamical decoupling pulse sequences. 
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Introduction. — Quantum systems are famously susceptible 
to interactions with their surrounding environments, a process 
which leads to a progressive loss of "quantumness" of these 
systems, via decoherence [ ]. When a system performs a 
quantum information processing (QIP) task this loss of quan- 
tumness is equivalent to the accumulation of computational 
errors, which leads to the eventual loss of any quantum advan- 
tage in information processing. Robust large-scale quantum 
information processing therefore requires that decoherence — 
or any otherwise undesired evolution — of a quantum state 
be minimized to the largest extent possible by clever sys- 
tem choice and engineering. One may then hope to apply 
the powerful techniques of fault tolerant quantum error cor- 
rection (FT-QEC) [2]. However, FT-QEC imposes signifi- 
cant resource requirements, in particular rapidly growing spa- 
tial and temporal overhead, together with demanding gate 
and memory error rates which must remain below a certain 
threshold (e.g., Refs. [i]). This motivates the search for al- 
ternative strategies which can slow down decoherence and 
"keep quantumness alive." Dynamical decoupling (DD) is a 
form of quantum error suppression that modifies the system- 
environment interaction so that its overall effects are very 
nearly self-canceling, thereby decoupling the system evolu- 
tion from that of the noise-inducing environment [ ]. DD has 
primarily been studied as a specialized control technique for 
quantum memory (i.e., arbitrary state preservation) [5-12], 
as convincingly demonstrated by a number of recent exper- 
iments in QIP platforms as diverse as electron-nuclear sys- 
tems [ij, i-f], photonics qubits [ ], and trapped ions [ ]. 
However, the holy grail of QIP is not just to store states ro- 
bustly, but rather to perform universal computation robustly 
[1/]. Fortunately, there are abstract results showing that DD 
is in principle compatible with computation, essentially by de- 
signing DD operations that commute with the computational 
operations [ ]. Additionally, recent theoretical results indi- 
cate that high fidelity "dynamically error-corrected gates" can 
be designed, using methods inspired by DD [ ], and that DD 
can be merged with FT-QEC to reduce resource overhead, or 
even improve gate error rates to below threshold [20]. DD 
can also be used to improve the fidelity of adiabatic quantum 
computation [ ]. Experimentally, DD has been successfully 
combined with QEC in nuclear spin systems to demonstrate 



robust quantum memory [^'']. In principle, then, it appears 
that DD is a suitable control technique for overcoming deco- 
herence and improving gate fidelity. However, a very practical 
question still remains: what are the conditions under which 
DD can be used to perform universal quantum computation 
with a given fidelity? Recent rigorous bounds devised for 
the popular "periodic DD" (PDD) protocol suggest that it is 
severely limited in this regard [ ]. Here we demonstrate, us- 
ing numerical simulations of a logical qubit coupled to a small 
bath, that recursively constructed, concatenated DD (CDD) 
pulse sequences [ ], can be used to endow a universal set of 
quantum logic gates with remarkably high fidelities. Though 
our numerical results illustrate the effectiveness of CDD in a 
model inspired by quantum dot implementations of QIP [ - ], 
the framework we describe is in principle more generally ap- 
plicable, provided encoding and encoded gates can be imple- 
mented. We therefore expect that our results will contribute to 
establishing DD as an indispensable tool in scalable QIP. 

Dynamical decoupling. — The total Hamiltonian without 
DD is H — Hb + HsB, where Hb includes all bath- 
only terms and Hsb includes all terms acting non-trivially 
on the system. To suppress error, DD allows the joint 
evolution to proceed under H for some time before ap- 
plying a control pulse Pj to the system alone [generated 
by a time-dependent system-only Hamiltonian Hg{t) which 
is added to H], designed to refocus the evolution toward 
the error-free ideal, continually repeating this process un- 
til some total evolution has been completed: DD[C/(to)] = 

Pjvf/(ro)---P2C/(ro)Pi?7(ro) = C/(iVro), where U{to) - 
Uq{tq)B{tq) represents the joint system-bath unitary evolu- 
tion generated by H, for a duration of length tq (the pulse 
interval), decomposed so that Uo{to) determines the ideal, de- 
sired system-only eiTor-free evolution, and B{to) is a unitary 
error operator acting jointly on the system and bath. Here 
and below we use a tilde to denote evolution in the presence 
of DD pulses. DD schemes with non-uniform pulse inter- 
vals have also been considered and shown to be very pow- 
erful for quantum memory purposes [9-12, 14, 16], but since 
it is unclear how to use them for computation involving mul- 
tiple qubits, we shall limit our discussion to uniform-interval 
schemes. For now, but not in our simulations presented be- 
low, we assume for simplicity of presentation that the pulses 
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{Pj} are sufficiently fast as to not contribute to the total time 
of the evolution. The simplest example is quantum mem- 
ory, where Uo{to) = Ig is the identity operation, and B{tq) 
represents the deviation from the ideal dynamics caused by 
the presence of a bath. In this case, our goal is to choose 
pulses so that DD[?7(ro)] = Is ^ B, where B is an arbi- 
trary pure-bath operator Uniform-interval DD schemes differ 
in precisely how the pulses {Pj} are chosen, with the only 
common constraint that the following basic "decoupling con- 
dition" (vanishing average Hamiltonian, i.e., vanishing first 
order term in the Magnus series of the joint system-bath evo- 
lution) is met [4]: X^a^Q^SB-^a = 0. To be concrete, we 
will suppose that the pulses Pa S {/, X, Y, Z} are Pauli 
operators. CDD generates pulse sequences by recursively 
building on a base sequence (motivated be- 

low), where [•] denotes either free evolution or the insertion 
of gate operations between pulses. The sequence is initialized 
as CDDo[C/(ro)] = [/(tq) - C/o(to)B(to) = C/o(to), and 
higher levels are generated via the rule CDD„_|_i[[/(to)] = 

Z[UniTn)]X[UM]Z[UniTn)]X[Un{Tn)] = C/„+i(t„+i), 

where = 4"to. Note that in contrast to previous work 
on CDD [6, 7, 20], we are allowing for the possibility of 
some non-trivial information processing operation Uo{to), as 
this will be required in our discussion of universal compu- 
tation below. The choice of the base sequence is motivated 
by the observation that it satisfies the "decoupling condi- 
tion", in the quantum memory setting Uq{to) = Is, under 
the dominant "1 -local" system-bath coupling term H^gg — 
E.e{...y,.}E,'^? ® Bf, where ^ X, a] ^ F, and 
(t| = Z denote the Pauli matrices acting on system qubit 
and {B"} are arbitrary bath operators. (The next order "2- 

local" coupling would have terms such as cr"cr^ ® ^i?f ' 
Similarly, the most common pulse sequence used thus far in 
DD experiments (e.g., Refs. [ , ]) is PDD, which generates 
pulse sequences by periodically repeating the base sequence 
PDDfe[t/(To)] = (PDDi[C/(to)])'= = 
J7fc(4fcro), where PDD (to)] = CDDi[C/(ro)]. Rigorous 
noise reduction bounds are known for both PDD and CDD 
in the quantum memory setting, and show that CDD is a 
much more effective strategy than PDD, provided (H-ffsH + 
ll^sslD'H] is sufficiently small, where the norm is the largest 
eigenvalue [ ]. It is convenient to characterize the leading- 
order DD behavior in terms of HifsH and ||i?ss||^ as these 
parameters capture the strength or overall rate of the internal 
bath and system-bath dynamics, respectively. If H-ffssll ^ 
then the system-bath coupling is a dominant source 
of error. In this case, DD should produce significant fidelity 
gains as it removes the dominant error source. On the other 
hand, if ||i?5B|| < H^^bH, then the system-bath coupling 
induces relatively slow dynamics, while the environment it- 
self has fast internal dynamics. In this case, suppressing the 
system-bath coupling will have less of an effect on the overall 
dynamics, so it may be considered a worst case scenario when 
assessing DD performance. 

Hi^h fidelity universal quantum gates using CDD. — Our 
main goal in this work is to demonstrate that we can gen- 
erate a universal set of logic gates which is highly robust 



in the presence of a decohering environment. As a model 
system we consider electron spin qubits in semiconductor 
quantum dots [24], which we study numerically via full- 
quantum-state (sometimes called "numerically exact") simu- 
lations over a wide range of system-bath coupling parame- 
ters. In such systems the dominant bath is provided by the 
nuclear spins p], and the interaction between system and en- 
vironment is described by a Heisenberg exchange Hamilto- 
nian with exponentially decreasing strength as a function of 
distance dij between system qubit j and bath qubit i. Thus, 
we let = J^iCrf /'^'^'^ in the system-bath Hamilto- 
nian H'^g^, so that ||i75B|| J. We model the interaction 
between the bath nuclear spin qubits as dipole-dipole cou- 
pling, i.e., Hb - /?E.<, (^f^^l + - 2af (Tj) /d% so 
that \\Hb\\ oc P, where now dij is the distance between bath 
qubits i and j. In our simulations we pick the parameters J, 
j3, and dij, as well as the pulse interval tq and the pulse width, 
to include a range of interest for GaAs and Si quantum dots 
[ , ]. The HsB and Hb Hamiltonians are on during the 
entire pulse sequence execution, while Hs{t) pulses appro- 
priately between dynamical decoupling {Hs = Hod) and 
computational operations {Hs — Hq). 

Universal quantum computation requires that only a dis- 
crete set of universal gates be implemented; a particularly 
simple choice are the Hadamard, 7r/8, and controUed-phase 
gates [17]. The first two are single-qubit gates, and the third 
is a two-qubit gate which can be used to generate entangle- 
ment. A conundrum immediately presents itself when trying 
to combine computation with DD: how to make sure that the 
DD pulses do not cancel the (system Hamiltonian implement- 
ing the) gates? One solution is to use an encoding so that 
the DD operations commute with the logical gate operations 
[ ! S, 21, 27]. To this end we use logical qubits encoded into a 
four-qubit decoherence-free subspace (DFS). The logical ba- 
sis states are the two orthonormal total spin-zero states of four 
spin-1/2 particles, first described as a DFS in Ref. [_b]. We 
stress that our system-bath interaction does not exhibit any 
symmetries so that there is no naturally occuring DFS which 
can be used to store protected quantum information; instead, 
our encoding choice is motivated by the fact that in this set- 
ting, a universal set of encoded computational operations can 
be generated by controllable Heisenberg exchange Hamilto- 
nians between the system qubits {— Hq), as first described 
in [ ], and these commute with the global Pauli operations 
{X,Z} — 1X1X2X3X4^, Z1Z2Z3Z4} used as decoupling 
pulses. To generate these pulses Hdd is modeled as a con- 
trollable uniform magnetic field. However, we emphasize that 
these choices are by no means unique. Any choice of DD 
pulses {Pj} such that [Hq, Pj] = Vj will suffice [' ], in- 
cluding, e.g., the stabilizer quantum error correcting codes rel- 
evant in the theory of FT-QEC used as DD pulses, and the nor- 
malizers of these codes used to generate computational gates 

Faced with several options for combining DD and compu- 
tational operations [ , ' ], we chose the following "decouple 
while compute" strategy. In this strategy we alternate between 
applying computational and DD operations, thus spreading a 
computational gate over the entire CDD pulse sequence. We 
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FIG. 1 : (color online) Fidelity of a universal set of encoded gates un- 
der CDD. The coupling strengths and bath dynamics are determined 
by the parameters J = lOkHz and /3 = IMHz, respectively (we 
work in units of = 1). Pulse intervals are fixed at tq = Ins, while 
pulse widths are given hy 5 = Q, 5 = Ips, and 5 = Ins, correspond- 
ing, from bottom to top, to the blue, green (absent for controlled- 
phase), and magenta lines, respectively. The red dashed line shows 
the unprotected evolution over a time period T — 4"tq. Notice that 
the logjQ(l — F) ranges change between plots. Also, n — corre- 
sponds to free evolution for a duration to = Ins, whence the n = 
point starts at a relatively high fidelity. Results depend only slightly 
on the choice of initial system state. 



do this by applying the N root of the gate A'^ times during 
a CDD pulse sequence involving N = 4" pulses. Thus, if 
the ideal computational gate is G(T) = e^*^<^^, we imple- 
ment it by applying U{to) = (,-i{Hg+Hsb+Hb)to between 
each of the N pulses, where tq = T/N . As T increases with 
concatenation level n, the exchange couplings in Hq are pro- 
portionately decreased so that the || i/cHT product remains 
constant; in all simulations the pulse interval tq is held fixed 
at Ins, the time-scale for exchange operations in semiconduc- 
tor quantum dots [ ]. This "decouple while compute" strat- 
egy is precisely the formulation presented in the expressions 
for CDD„+i [J7(ro)] and VDT) k\U [tq]] above, provided we 
identify [/o(To)i?(To) there with U (tq) here, and ?7o(to) there 
with G{tq) here. Other strategies are certainly also conceiv- 
able, e.g., a "decouple then compute" strategy wherein Uq{tq) 
is simply the identity operation, and the gate is implemented 
at the end of the pulse sequence. While the latter strategy was 
shown to be capable of reducing the resource requirements of 
FT-QEC [ ], we found in our simulations that we obtain a 
higher fidelity when we use "decouple while compute", be- 
cause then time is not wasted on free evolution during the in- 
tervals between pulses. 

We now present our simulation results (details of the nu- 
merical procedure are given in [ ]). The worst-case scenario 
of J < is shown in Fig. 1 , where we plot log]^g(l — F) vs. 
concatenation level for each of the universal gates, with the 
fidelity defined as = -\/| p |, where p is the mixed 
output system state (obtained from the joint system-bath evo- 
lution after partial trace over the bath) and is the desired 
system state. In each of these plots, the red dashed line rep- 



resents undecoupled free evolution for increasing total time, 
given by T = 4'Vo. As the evolution time increases, error ac- 
cumulates and fidelity correspondingly worsens, while CDD„ 
combats this effect with each successive level of concatena- 
tion. To contrast, the blue line in these graphs shows CDD„ 
with ideal, zero-width DD pulses, so that realistic, finite-width 
DD pulses lie somewhere between the blue and red lines, as 
shown. In each of the plots in Fig. 1 CDD achieves impressive 
results, even when pulse widths are as long as the intervals, 
that is, when 5 = tq = Ins as depicted with the magenta lines, 
CDD still manages more than five orders of magnitude im- 
provement in fidelity over free evolution. As the pulse width 
5 narrows relative to the pulse interval tq, that is, as the DD 
pulses becomes faster, fidelity improvement grows to between 
ten and twenty orders of magnitude over free evolution. 

The results for the encoded 7r/8 and Hadamard gates are 
similar, which is not surprising given that they require, re- 
spectively, one and two elementary Heisenberg exchange op- 
erations to be implemented [ , ]. The fidelity of the 
controlled-phase gate is several orders of magnitude lower, 
which is due to the fact that it involves a much longer sequence 
of 42 elementary Heisenberg operations [ ]. Finally, while 
the quantum memory results are comparable to those of the 
Hadamard and 7r/8-gates, we attribute the reduction in mem- 
ory fidelity at the highest concatenation level to the absence 
of Hg during the intervals between pulses. Indeed, having 
the system Hamiltonian "on" during the pulse intervals has a 
beneficial effect, as it effectively reduces the strength of the 
bath and system-bath Hamiltonians. The overall conclusion 
from Fig. 1 is rather encouraging: it appears to be possible to 
implement a universal set of quantum logic gates with a high 
fidelity in the presence of coupling to a spin bath. 

The results in Fig. 1 are for specific coupling parameters 
chosen deliberately to represent a worst-case scenario for DD, 
in that J < /3. As we next demonstrate, the conclusions are 
robust: CDD remains effective over a broad range of bath dy- 
namics and system-bath coupling strengths. Figure 2 shows 
the resilience of CDD to widely varying environments by dis- 
playing constant fidelity contours in ( Jtq, /3ro) space, at fixed 
concatenation level and pulse width, as indicated. Note that 
fixing n and 5 renders the total evolution time constant, so 
that fidelity becomes strictly a function of the dimensionless 
coupling parameters (Jro,/3ro). These plots show a strong 
fidelity dependence on Jtq, and a very weak dependence on 
/3ro, except in the quantum memory case. 

More generally, our results show that CDD is effective over 
a broad range of coupling parameters, including the funda- 
mentally different "good" (J > /3) and "bad" (J < P) 
regimes. This conclusion is further bolstered by our complete 
gate fidelity simulations [ ], where the /3 and J parameters 
each vary over the range from IHz to IMHz. In these simu- 
lations the non-memory gate fidelities improve monotonically 
as a function of concatenation level for all values of J and j3. 
Taken in their totality, our simulation results indicate that uni- 
versal quantum computation can be combined with CDD to 
achieve very high fidelities. 

Discussion. — The high gate fidelities we have reported here 
suggest that it is advantageous to incorporate CDD as a first 
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FIG. 2: (color online) Constant fidelity contours for the system described in the previous figure, at fixed concatenation level n = 5 and pulse 
width 5 — Ins. Notice that the fidelity contours are strongly dependent on Jtq. 



layer of defense against decoherence, in a more complete FT- 
QEC scheme. In this regard our "decouple while compute" 
study complements the "decouple then compute" strategy for 
which it has already been shown that incorporating DD, and 
in particular CDD, into FT-QEC can lead to substantial im- 
provements [20]. It will be interesting to determine which 
strategy leads to better performance overall. Undoubtedly, op- 
timal control methods will offer important additional perfor- 
mance improvements in the implementation of encoded logic 
gates [ "]■ We look forward to experimental tests of CDD- 



protected quantum memory and logic gates. 
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Appendix A: The DFS code 

We describe the four-qubit DFS code, first proposed in the context of providing immunity against collective decoherence 
processes [ ]. Let S and ms denote the quantum numbers associated with total spin and its projection, and let the singlet and 
triplet states of two electrons i,j he denoted as 
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Here denotes a normalized basis state with the first (second) electron in the spin up (down) state, etc. Then a single en- 

coded DFS qubit is formed by the two singlets of four spins, i.e., the two states with zero total spin — | ^" -|- -|- S/^ |, 
where Si is the Pauli spin vector-operator of electron i. These states are formed by combining two singlets of two pairs of spins 
(|0l)), or triplets of two pairs of spins (|1l)), with appropriate Clebsch-Gordan coefficients: 

|0l) = \s)ab(S\s)cd 

= \ (i*(nn)> + i*(;nt)> - i*(mt)> - i*(;m)» 

= -L(2|M^(tm)> + 2|v&(;;n)> - i*(mt)> - i*(;m)> - i*(nn)> - i*(;nt)». 

The details of the sequences of exchange interactions needed to implement the Hadamard, 7r/8, and controlled-phase gates 
over this DFS code are too long to give here. However, they have been well documented in the original works [-)_-.i-r], and 
perhaps most concisely in the more recent Ref. [35]. 



Appendix B: Details of the numerical procedure 

We performed numerically exact simulations since we required extremely high precision fidelity results, inaccessible via 
approximation techniques [36, 37] capable of handling much larger system and bath sizes. Such methods are required for 
long time scale simulations, which is not our case. Moreover, large-scale simulations [ ] have confirmed earlier small-bath 
simulations of CDD pulse sequences [38, 39]. Our numerically exact simulations ran at 100-200 digits of numerical precision, 
for more than 90 hours on a computer with a dual core intel processor (2GHz + 2G RAM). In our simulations, each logical qubit 
was encoded using four physical qubits. We tested for dependence on the initial encoded system state, and found a variation in 
output fidelity of less than an order of magnitude. Hence all our reported results are for the logical-one irutial state. As a further 
test we considered several simple bath geometries: linear, circular, and polygonal. Our results did not depend appreciably on 
this geometry. We took the initial bath state as the uniform superposition, zero temperature state -j= N)' where B = 2^" 
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is the number of available pure bath states when Nn bath qubits are present. We checked that finite temperature has only a small 
quantitative effect on our reported fidelities. In order to keep our numerically exact simulations feasible we had to restrict the 
total number of qubits to ten, and hence the size of the bath to only two physical qubits when considering two encoded 4-qubit 
DFS qubits. For consistency we also used two bath qubits in our single encoded qubit simulations. However, we adjusted the 
strength of the coupling to the bath to account for this, using a scaling relation we found by testing bath sizes in the range 
2 < Nn < 5. While we are fully aware of the potential problems associated with using a small bath at long time scales, we 
expect our small bath simulations to be reliable, as we are concerned only with short timescales, relative to the magnitudes of J 
and /3, so that possible coherence recurrences due to non-Markovian effects are well out of reach. 



Figures 3-8 present our complete gate fidelity results for the 4-qubit DFS code subject to free evolution or CDD, with the j3 
and J parameters each varying over the range from IHz to IMHz. The pulse interval tq is fixed at Ins, and the pulse width 
assumes the values 0, Ips, and Ins. In all cases CDD leads to a fidelity improvement of many orders of magnitude relative to 
free evolution. The plots surrounded by dashed line boxes are the ones shown in Fig. 1 in the main text. 



We picked the ranges of our j3 and J parameters shown in Figs. 3-8 to represent realistic quantum dot systems - see Table 1. 
The J parameters in this table were estimated directly from the contact hyperfine interaction between a quantum dot electron 
and nuclear spins. For GaAs the hyperfine constant is yl = 90/zeV [40]. For polarized nuclei, the interaction strength is 
reduced by the number of nuclei N^, typically in the range 10^ to 10^ per quantum dot, giving J — A/N^- For unpolarized 
nuclei, the interaction strength is reduced by \/Nn [ ]. For Si (silicon), the contact hyperfine strength A is computed from 
A = ^ iJ,29SiHB'''i'r], where nj is the nuclear or electron (Bohr) magnetic moment, r is the fraction of non-zero spin ^^Si 
nuclei, v is the Si nuclear number density, and rj is the concentration of the electron Bloch wavefunction near the Si nuclei 
[41]. For natural abundance Si we obtained A = 60neV. For polarized nuclei, the interaction strength is reduced by the number 
Nn ~ 6 X 10^ of ^^Si in a quantum dot. For unpolarized nuclei the interaction strength is again reduced by ^/Nn- 

For both GaAs and Si the energy scale /3 of the bath is taken to be 1/T2„, the inverse of the nuclear dephasing time. For GaAs, 
the nuclear dephasing time has been estimated at = 100/is [42, 43], which is the precession time of a nucleus in the dipole 
magnetic field generated by neighboring nuclei. For Si, the nuclear dephasing time has been measured to be T2n = 5.6ms [ ], 
again corresponding to the evolution time due to the dipole-dipole interaction between a ^^Si nucleus and neighboring ^^Si 
nuclei. 



Appendix C: Complete simulation results 



Appendix D: GaAs and Si parameters 



Table 1 



system 



J /3 



unpolarized GaAs 
polarized GaAs 
unpolarized Si 
polarized Si 



lOOMHz lOkHz 
IMHz lOkHz 
400kHz 180Hz 
1.5kHz 180Hz 
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FIG. 3: CDD quantum memory in the 4-qubit DFS encoding with pulse interval tq — 10^ s, and pulse widths S — Os, S = 10^^ s, 
and S — 10~®s corresponding to the blue, green, and magenta lines, respectively. The dashed-red represents free evolution without DD for 
increasing total times T — 4"tq. The plots are labeled by the strength of the coupling parameters (logj^Q J, logj^p Z^)- Hz. Thus J increases 
in strength in multiples of 10 from left to right, and /3 similarly increases from top to bottom. Notice that the logjo(l — Fidelity) ranges change 
between plots. 




FIG. 5: Top: Turning point of the memory fidelity shown in Figure 3, i.e., the maximum concatenation level before fidelity decreases, as a 
function of log^o J and log^g in Hz. Bottom: Same for the vr/S gate fidelity shown in Figure 4. 
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Concatenation Level (n) 

FIG. 6: Same as Figure 3, for the Hadamard gate. 




FIG. 7: Turning point of the Hadamard gate fidelity shown in Figure 6, i.e., the maximum concatenation level before fidelity decreases, as a 
function of logjQ J and logig /3, in Hz. Note that slightly better performance is obtained for the n/8 gate. This is due to the fact that it takes 
only a single elementary Heisenberg exchange operation to implement, while the Hadamard gate takes two such operations, i.e., double the 
time. 



